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SEMI-ANNUAL STATUS REPORT 

NASA Grant NGR 05-002-100 

Space Radiation Laboratory (SRL) 
California Institute of Technology 

1 April 1904 - 31 March 1905 


This report covers the research activities In Cosmic Rays, Gamma Rays, and 
Astrophyslcal Plasmas supported under NASA Grant NGR 05-002-160. The report is 
divided into sections which describe the activities, followed by a hi holography. 

This group's research program Is directed toward the investigation of the astro- 
physical aspects of cosmic rays and gamma rays and of the radiation and elec- 
tromagnetic Qeld environment of the Earth and other planets. We carry out these 
investigations by means of energetic particle and photon detector systems flown on 
spacecraft and balloons. 


1. Cosmic Rays and Astrophysical Plasmas 

This research program is directed toward the investigation of galactic, solar, 
Interplanetary, and planetary energetic particles and plasmas. The emphasis is on 
precision measurements with high resolution in charge, mass, and energy. The main 
efforts of this group, which are supported partially or fully by this grant, have been 
directed toward the following two categories of experiments. 


1.1. Activities in Support of or in Preparation for Spacecraft Experiments 

These activities generally embrace prototypes of experiments on existing or 
future NASA spacecraft or they complement and/or support such observations. 


I. 1 .1. The High Energy Isotope Spectrometer Telescope (HEIST) 

HEIST is a large area (0.25 m^sr) balloon-borne isotope spectrometer designed 
to make high-resolution measurements of isotopes in the element range from neon 
to nickel (10 £ Z ^ 20) at energies of about 2 GeWnucleon. The instrument consists 
of a stack o'. 12 Nal(Tl) scintillators, two Cerenkov counters (Cl and C2), and two 
plastic scintillators (SI and S2) as illustrated in Figure 1. Each of the 2-cm thick Nad 
disks is viewed by six 1.5-inch photomultipliers whose combined outputs measure the 
energy deposition in that layer. In addition, the six outputs from each disk are com- 
pared to determine the position at which incident nuclei traverse each layer to an 
accuracy of ~2 mm. The Cerenkov counters, which measure particle velocity, are 
each viewed by twelve 5-Lnch photomultipliers using light integration boxes. This 
experiment is a collaborative effort with the Danish Space Research Institute. 
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Figure 1 

HEIST determines the mass of Individual nuclei by measuring both the change in. 
the Lorentz factor (A7) that results from traversing the Nal stack, and the energy 
loss (AE) in the stack. Since the total energy of an isotope is given by E = 7M, the 
mass M can be determined by M = AE/A7. The instrument is designed to achieve a 
typical mass resolution of 0.2 amu. 

On May 14, 1964, the HEIST instrument was launched from Palestine, Texas on 
its first bailoon flight. It reached altitude successfully and floated over Texas for 
approximately 30 hours at. a typical atmospheric depth of '-'5.5 g/cm . and ai 
geomagnetic cutoffs ranging from ~4J5 to 5.5 GV. During the flight more than 4x10 s 
events were recorded, including ~10 £) events with charge Z i 8 and kinetic energy E 
a 1.4 GeV/nucieon. Although only a subset of these are suitable for isotope analysis, 
the others are being used for in-flight mapping and stability checks. 

The various HEIST subsystems operated near-perfectiy during the flight. 
Preliminary analysis of the flight data has shown that ail of the 108 photomultiplier 
tubes (PMTs) operated as expected, and that their outputs were stable to ^ 1% over 
the course of the flight. The level of stability was a direct result of the successful 
operation of the automatic evaporative cooling system, which maintained the 



internal temperature of the gondola to 25 ± 0.5°C, thereby minimizing the need for 
tomperature-dependent corrections to the gains of the detector electronics, or the 
Nal light output. 

During the past six months significant progress was made In developing 
algorithms for the determination of the position and energy loss of Ions traversing 
the Nal stack. These algorithms depend on detailed maps ml the response of each 
PMT viewing the Nai; such maps are being derived from the & °Mn Bevalac calibration 
of HEIST. 



Figure 2 

Using these maps in association with appropriate interpolation algorithms we 
have now achieved an rms energy-loss resolution of ~ t.2% in the first Nal layer, 
averaged over the entire area of the disk as shown in Figure 2. This value is only 
slightly worse than the value achieved (~ 1,0%) in limited regions of the center of the 
disk. To determine position in the Nal disks we have developed algorithms based on 
various ratios of the response of the six PMTs viewing each disk. We are now 
achieving a typical resolution of ~ 2 mm for &t) Mn ions incident over ~ 90% of the 
diameter of the first Nal disk, as illustrated in Figure 3. 
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Figure 3 

The position and energy resolution currently being achieved is consistent with 
the requirements for isotope resolution in HEIST. We are now In the process of 
mapping the Nal layers deeper in the stack so that the information from the various 
layers can be combined to derive the energy loss of each event, identify nuclear 
interactions, and determine the trajectory of individual events. 

Our colleagues at the Danish Space Research Institute (DSRI) have also received 
copies of the flight data tapes. They are concentrating their efforts on mapping the 
response of the two Cerenkov counters, and calibrating the methods for velocity 
determination. 

The following recent talks and papers have resulted from this work; 

•"Measurements of Cosmic Ray Isotopes from a Space Station Platform,"R. A. 
Mewaldt, Workshop on Cosmic Rays and High Energy Gamma Ray Experiments 
for the Space Station Era", Paper presented at the Workshop on Cosmic Rays 
and High Energy, Louisiana State University, October, 1984 (1984). 

• "Initial Results from the Caltech/DSRI Balloon-Borne Isotope Experiment." S. 

M, Schindler et al . Proc. 19th International Cosmic Ray Conference, La Jolla. 
California (1985). 

A summary of an Invited paper that was presented at the LSU Workshop on 
Space Station Experiments appears below. 
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MEASUREMENTS OF COSMIC RAY ISOTOPES FROM A SPACE PLATFORM 

The Isotopic composition of galactic cosmic rays contains a record of the 
nuclear history of a sample of matter from other regions of the galaxy, 
Including Its synthesis In stars, and its subsequent nuclear interactions with 
the interstellar gas, Although it Is only recently that high resolution 
observations of cosmic ray isotopes became experimentally possible, they have 
already altered our views of both cosmic ray origin and propagation, In spite of 
thl3 progress, the rich potential of these studies remains largely unrealized 
because hlgh-resotution observations to date have been statistically limited - 
only the most abundant species have so far been measured. Much of the Lmtial 
progress In this field was made at low energies (-30 to 300 MoWnucleon) with 
small solid-state spectrometers carried on spacecraft outside the 
magnetosphere, Much larger balloon-borne instruments have now been 
developed to resolve Isotopes at energies from -0.5 to -2 GeV/ nucleon. 
Recently, consideration has been given to the feasibility of a super-conducting 
magnet facility for the space station, that, among other objectives, might 
measure cosmic ray Isotopes up to energies as high as -100 GeV/nucleon. In 
order to capitalize on these new capabilities and make statistically accurate 
measurements of rare Isotopic species, It will be necessary to expose large 
Instruments outside the atmosphere for periods of a year or more. Such 
experiments are therefore obvious candidates for a "space platform", either 
associated witn the space station itself, or a "freo-flyer" launched into a 
hlgher-lnclinatton orbit. In this paper various possibilities for measuring 
cosmic ray isotopes with 3=:Z^3Q from an Earth-orbiting space-platform are 
■•onsldered. It is found that space platform experiments could provide a factor 
of — 100 improvement in yield over currently approved future experiments. To 
address the objectives of cosmic ray isotope studies most effectively, 
measurements should be made by several complementary Instruments spanning 
a range of energies from <0,1 to -100 GeV/nucleon. 

Objectives and Requirements of Cosmic Ray Isotope Studies: Some of the 
principal objectives of cosmic ray isotope studies Include: (l) the study of 
nucleosynethesis in other regions of the galaxy, (2) measurement of the time- 
delay between nucleosynthesis and cosmic ray acceleration, and (3) studies of 
the cosmic ray storage time and pathiength distribution In the galaxy. To 
address these objectives will require precision measurements of a variety of 
isotopes with l^Zi30, including 1) rare, neutron-rich Lsotopes of C, 0, Ne, Mg, 
Si, S, Ca, Fe, and Ni; (2) electron-capture Isotopes such as a7 Co; and (3) 
radioactive clocks such as 10 Be, 20 Al, and 3a Cl. These measurements must have 
sufficient mass resolution (^0,25 amu) and collecting power (e.g,, a 10° SI 
events), preferably in several different energy ranges. For comparison, 
Instruments flown to date have collected a total of fewer than 1000 Si (Z = 14) 
events with (T m ^ 0.3 amu. However, the required measurements are now 
experimentally possible if present-day instruments could be exposed In space 
for sufficient length of time. 
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Figure 4 

Element and Energy Coverage: Figure 4 summarizes the element and energy 
ranges covered by representative present-day isotope spectrometers. At the 
lowest energies the energy loss - total energy (AE-E‘) technique has been very 
successful, but is best suited for spacecraft outside the magnetosphere. 
Except for the magnet experiment, the remaining methods are all based on 
Cerenkov counters, including the Cerenkov - range (C-R), Cerenkov - total 
energy (C-E), and Cerenkov - energy loss - Cerenkov (C-E-C) techniques. Note 
that none of these methods can presently extend above GeV/nucleon. 



- 9- 


or r-'cjc:,* 


r. 


'S 


OUAL/TV 


FUTURE ISOTOPE SPECTROMETERS 



Figure 5 

Figure 5 shows estimates of the coverage for future Isotope spectrometers. 
Each of the various Cerenkov techniques can to some extent be '’tuned 1 ' in 
energy by selecting the Index of refractLon of the Cerenkov radiator(s), which 
allows for some flexibility In .he energy and charge ranges covered by a 
particular Instrument (note that a single instrument configuration would not 
necessarily cover the entire energy range Indicated). Also shown in Figure 5 
are two possible designs for super-conducting magnet experiments that 
combine a measurement of the particle's rigidity (provided by the magnet and a 
trajectory measuring device) with a measurement of the particle’s velocity (or 
momentum/nucleon), perhaps from a Cerenkov counter. The upper limits to the 
energy range Indicated for the "shuttle-compatible" and "space station" 
versions were determined solely from the assumed rigidity resolution of the 
magnet; assuming that velocity measuring devices of sufficient accuracy could 
be developed (undoubtedly a considerable experimental challenge, since to 
resolve Fe isotopes with <7 m = 0.25 amu requires that the momentum/nucleon 
be measured to ~0.4 % accuracy or better), It Ls Likely that any given velocity 
measuring device would function over only a limited portion of the energy 
range accessible to the magnet. Although a superconducting magnet facility 
represents a considerably more ambitious undertaking than the "Pi-class" 
experiments in, Figures 4 and 5, it offers the possibility of extending cosmic ray 
isotope measurements into the energy region from ~10 to ~100 GeW nucleon. 
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The Collecting Power of Earth-Orbiting Spec trocne tors: Tile principal objoctlvo 
in placing a cosmic ray Isotope spectrometer in orbit on a space platform is to 
Improve its yield so that the statistical uncertainty of measurements of rare 
Isotopic species can be minimized. By far the most important gain by a space 
platform oxperimrnt Is in the exposure lima. None of the existing balloon 
isotope expertmr' .s have been exposed for more than a day or two; this time 
should bo more than a factor of 100 greater on a space platform. The primary 
disadvantage of an Earth-orbiting platform Is the affect of the geomagnetic 
cutoff. 
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Figure 6 

Figure 9 shows the time-averaged spectra that would be obsrrved by an Earth- 
orbiting cosmic ray experiment for several orbital inclinations. Some of the 
orbits shown have special significance. The 28.5 deg orbit corresponds to the 
latitude of Cape Kennedy, and it will undoubtedly be the orbit of the first space 
station. An orbital Inclination of 57 deg Is the maximum obtainable with a Cape 
Kennedy shuttle launch. By launching from Vandenberg APB on the west coa3t, 
a polar orbit (\ = 90 deg) can be obtained. Note that at energies >7 
GeV/nucleon the geomagnetic cutoff has no effect. However, using present 
techniques (including the magnet approach), cosmic ray isotopes are most 
easily measured at somewhat lower energies, as shown In Figures 4 and 5, In 
the energy range from ~0.5 to ~2 GeV/nucleon it is clearly Important to have a 
high-inclination orbit (A^45 deg), or much of the benefit of the increased 
exposure time gained by putting an experiment into orbit will be lost. 

Yield Estimates: Figure 7 shows how the collecting power of three possible space 
platform experiments depends on the inclination of the orbit, and also 
compares these to various (lower- energy) spectrometers on interplanetary- 
spacecraft that have either fiown, or are approved and/or proposed. Note that 
the C-R, C-E, and C*E'C methods ali require a relatively high- inclination orbit to 
be effective. An advantage of the magnet approach is that it can function at 
higher energies, and thus is compatible with a 28.5 deg orbit such as the first 
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space station will have, although In “ oh an orbit it could not contribute 
measurements below ~2 GeV/nucloon (see Figure 0) Note that each of these 
three experiments In a suitable orbit (and also an Interplanetary Explorer) 
could easily achieve a factor of ~100 to 1000 Increase In yield ovor experiments 
that have flown to hate (Including ISEE-3 and all reported balloon experiments), 
and they could also Improve by a factor of -100 ovor presently approvod 
missions such as ISPM (now named Ulysses). 



Figure 7 

5 

As Figure 7 shows, it Is not unreasonable to expect yields of a few times 10 
SI events per year In several energy ranges if these or similar experiments 
could be flown in a suitable orbit, Indeed, ytelds as large as ~10° Si per year 
are possible if the instruments were scaled up further In size over the next few 
years, 



Figure B 

Figure 8 shows the corresponding yield of other isotopes with 3^5:30, for SI =< 3 
x 10 . Note that for most of the isotopes of interest, including the radioactive 
clocks, the expected yield ranges from ~10^ to -10“ per year, sufficient to 
make definitive measurements of even rare species. 

Sunmury In summary, there are now techniques available that can resolve 
cosmic ray Isotopes with Zi30 over a wide range In energy, from <0.1 
GeV/nucleon up to perhaps even 100 GeWnucleon. If flown on a space platform 
such Instruments could improve over present or approved experiments by at 
least two orders of magnitude In yield. In addressing the objectives of cosmic 
ray Isotope studies it will be essential to measure the composition In several 
different energy ranges. In addition to the time-dilation effects on radioactive 
clocks, there are other energy-dependent processes, including acceleration by 
supernova shock waves, fragmentation, solar modulation, and apparently the 
galactic confinement time, all of which will require broad energy coverage to 
understand. A super-conducting magnet facility on the space station would 
make a major contribution to this need by extending Isotope measurements to 
very high energies, thereby making accessible a new range of time-dilation 
factors for radioactive clocks, and providing a measure of the composition in an 
energy region where fragmentation corrections are greatly reduced so that the 
source abundances of a greater number of Isotopes could be accurately 
determined. However, since space station observations will be Limited to 
energies greater than —2 GeV/nucleon, well above those of the bulk of cosmic 
ray nuclei, It will be necessary to launch some experiments Into higher- 
IncUnatlon orbits and/or outside the magnetosphere if the required energy 
coverage is to be obtained. For this reason it important that we develop a 
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number of complementary experimental approaches to measuring cosmic ray 
Isotopes, because It Is essential that those studios span ad broad an onergy 
range as posnlblo Lf we oro to road accurately the record that cosmic ray 
Isotopes transmit. 


1.1.2. Low Energy Isotope Spectroscopy 

In August of 1904 we completed a successful calibration of solid-state 
detectors at the Lawrence Berkeley Laboratory Bovaiac, in which wo studied 
scvoral problems Important to isotope spoctromotors based on silicon solid- 
state detectors. In particular, wo investigated a "charge multiplication” 
problem that has plagued several thin (~30 to 150 p.) solid-state detectors flown 
on previous missions. The characteristic of this phenomenon is that the signal 
from a significant fraction of the portlclos that lose a great deal of energy in 
the dotoctor (e.g., stopping Fe tons) Is measured to be ~20 to 30>4 greater than 
expected, wtth greatly degraded resolution. At Berkeley wo tested a detector 
from ORTEC that had received special maunufacturing treatments designed to 
eliminate this problem. Preliminary analysis indicates at least an order of 
magnitude reduction in the fraction ol events exhibiting this problem. 

Other areas of investigation during this calibration were an extension of 
our mass spectroscopy capabilities up to Kr (Z=38) isotopes, further 
calibrations of pos'Uon-sensiUve solid-state detectors, and further work on 
determining cosmic ray fragmentation cross sections at low energies. These 
data are now being analyzed. 


1.2. Experiments on NASA Spacecraft 

The SR&T grant program of the Space Radiation Laboratory is 
strengthened by and contributes to the other programs described here. 
Activities related to these programs are primarily funded by mission- related 
contracts but grant funds are used to provide a general support base and the 
facilities which make these programs possible. 


1.2.1. An Electron/laotope Spectrometer (F3S) Launched on IMF-7 on 22 
September 1272 and on IMP-B on 26 October 1973 

This experiment Is designed to measure the energy spectra of electrons 
and positrons (0,16 to ~8 MeV). and the differential energy spectra of the nu- 
clear isotopes of hydrogen, helium, lithium, and beryllium (~2 to 50 
MeV/nucleon), In addition, It provides measurements of the fluxes of the Iso- 
topes of carbon, nitrogen, and oxygen from ~5 to ~15 MeV/nucleon. The meas- 
urements from this experiment support studies of the origin, propagation, and 
solar modulation of galactic cosmic rays; the acceleration and propagation of 
solar flare and interplanetary particles; and the origin and transport of ener- 
getic magnetospheric particles observed in the plasma sheet., adjacent to the 
magnetopause, and upstream of the bow shock. 
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The extensive EIS data set has been utilized In comprehensive studies of 
solar, Interplanetary, and magnetosphertc processes. Correlative studies have 
involved data from other IMP Investigations and from other spacecraft, as well 
as direct comparisons of EIS data from IMP-7 and IMP-0. In partlci ler, we are 
currently using IMP data to support a multidisciplinary workshop that Is study- 
ing cosmic ray modulation and the largo scale structure of the heliosphere, 

Our studies of IMP data have resulted in the following rocent papers; 

• "Microstructure of Magnetic Reconnection in Earth’s Magnetotail," J. W. 

Bleber ot al., J. Geophys, ffcs. 09, 0705 (1904). 

■ "The Isotopic Composition of the Anomalous Low-Energy Cosmic Rays," 

R.A. Mewaldt et al„ /Ip. /. 263, 450 (1904), 

• "Solar Cycle Variations of the Anomalous Cosmic Ray Component." R A. 

Mewaldt and E. C. Stone, Proc. 19th international Cosmic Ray Conference, 

La Jolla, California (1905), 

q 

• "Cosmic Ray He Measurements," R. A Mewaldt. Proc. 19t.h International 

Cosmic Ray Conference, La Jolla, California (1905). 


1 .3.2. An Interstellar Cosmic Ray and Planetary Magnetospheres Experiment for 
the Voyager Missions Launched In 1977. 

This experiment is conducted by this grou.p m collaboration with F. B. 
McDonald and J. H. Trainor (Goddard Space Flight Center), W. R. Webber (Univer- 
sity of New Hampshire), and J. R. Jokipti (University of Arizona), and has been 
designated the Cosmic Ray Subsystem (CRS) for the Voyager Missions. The ex- 
periment is designed to measure the energy spectra, elemental and (for lighter 
elements) isotopic composition, and streaming patterns of cosmic -ray nuclei 
from H to Fe over an energy range of 0,5 to 500 MeV/nucleon and the energy 
spectra of electrons with 3 - 100 MeV, These measurements will be of particular 
importance to studies of stellar nucleosynthesis, and of the origin, acceleration, 
and interstellar propagation of cosmic rays. Measurements of the energy spec- 
tra and composition of energetic particles trapped in the magnetospheres of 
the outer planets are used to study their origin and relationship to other physi- 
cal phenomena and parameters of those planets. Measurements of the intensi- 
ty and directional characteristics of solar and galactic energetic particles as a 
function of the heliocentric distance will be used for in sthu studies of the in- 
terplanetary medium and its boundary with the interstellar medium Measure- 
ments of solar energetic particles are crucial to understanding solar composi- 
tion and solar acceleration processes. 

The CRS flight units on both Voyager spacecraft have been operating suc- 
cessfully since the launches on August 20, 1977 and September 5, 1977. The 
CRS team participated in the Voyager 1 and 2 Jupiter encounter operations in 
March and July 1979, and in the Voyager 1 and 2 Saturn encounters in No- 
vember 1900 and August 1 98 1 . The Voyager data represent an immense and 
diverse data base, and a number of scientific problems are under analysis. 
These investigation topics range from the study of galactic particles to particle 
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acceleration phenomena In the Interplanetary medium to plasma /field ener- 
getic particle Interactions, to acceleration processes on the sun, to studies of 
elemental abundances of solar, planetary, interplanetary, and galactic energet- 
ic particles, and to studies of particle /field /satellite interactions In the magne- 
tospheres of Jupiter and Saturn. 

The following publications and papers for scientific meetings, based on 
Voyager data, were generated: 

• "Elemental Composition of Solar Energetic Particles," W.R. Cook et al ., 
Ap. J. 279, B27-B38 (1904), 

• "The Voyager Mission: Encounters with Saturn,” E.C. Stone, J. Geophys, 
Res. 60, 8369 (1983). 

• ,r Voyager Measurements of the Energy Spectr um. Charge Composition, 
and Long Term Temporal Variations of the Anomalous Components In 
1977-1982,” W. R. Webber and A. C. Cumming3, Proceedings of Solar Wind 
V Conference. NASA CP 2280, 427-434 (1983). 

• ‘ Temporal Variations of the Anomalous Oxygen Component," A. C. Cum- 
mings and tf. R. Webber, Proceedings of Solar Wind V Conference, NASA CP 
2280,435-440(1983). 

• “The Saturn System." E.C. Stone and T. Owen, Shh/m, 3 (1984). 

• "E'hdence that the Anomalous Component is Singly Ionized." A. C. Cum- 
mings et al., Ap. J. (Lett.) 287, L99 (1904). 


One of these studies is summarized here: 

Latitude Variation of Recurrent Protons in the Heliocentric Radial Range 

11 to 20 AU and Correlation with Solar Coronal Hole DynanAcs. 

During the decline to minimum sunspot activity of the last 11-year solar 
cycle (1973 to 1977) recurring enhancements of MeV-energy nucleon fluxes 
were found to be closely associated with plasma interaction regions corotating 
with the sun near the heliographic equator at distances of 1 to 8 AU from the 
sun. These "corotating interaction regions' 1 , or CIRs, develop at the interface 
between long-lived high speed solar wind streams and lower speed solar wind 
streams. Recurrent flux enhancements most likely result from the interplane- 
tary acceleration of keV-energy nuclei up to MeV energies by multiple interac- 
tions with the hydromagnetic shocks that often form on the leading and trailing 
edges of the CIRs beyond 1 AU. Elemental composition studies indicate that the 
keV-energy seed population is most likely the suprathermal tail of the solar 
wind, although solar ftxre accelerated nuclei which are also present can be 
further accelerated along with the interplanetary seed population. Recurrent 
flux enhancements have previously been used to indicate the presence of the 
CIRs. 

The sources of high speed solar wind streams are now known to be coronal 
holes, the longest lasting of which cover the poles of the sun during activity 
minimum. Equatorial extensions of the polar holes are often sources of high 
speed streams observed near the heliographic equator. Following the peak sun- 
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spot activity in late 1979 to mid 1980, long-lo.3tl.ng holes reappeared In the sun's 
corona, allowing Investigations of the association between coronal hole and 
solar wind stream dynamics and the Interplanetary acceleration of nuclei to 
MoV energies to be extended Lnto th.6 outer solar system using data from Ihe 
Voyager and Pioneer spacecraft, 

Spacecraft Observations of Recurrent Phenomena: This study is based on daily 
samples of the average counting rate of protons In the energy range ~0,5 to 17 
MeV from Low Energy Telescopes of the Voyager CRS and the counting rate of 
~0.B to 1.0 MeV protons from the GSFC-UNH Cosmic Ray Experiment aboard the 
Pioneer 11 spacecraft. Spacecraft distance from the sun (R). latitude (9), and 
longitude ($) In an inertial heliographic coordinate system are listed in Table 1 
for the start and end of the analysts period. Pioneer 1 1, midway between the 
Voyagers in radius, Is closer to Voyager 2 In longitude (<7° separation) but is 
closer to Voyager 1 in latitude (~6° to 9° separation). 



Voyager 2 

Pioneer 11 

Voyager 1 

R 

11.0- 14,2 AU | 

13.1 - 18.1 AU 

15.1 - 

19.8 AU 

0 

2.4°S - Q.8°S 1 

1 1.7°N - 14.3°N 

17. 1°N 

- 23.1°N 

4> 

I 327,1° - 348.4° ' 

329.6° - 342.6° 

315.6° 

- 325.7° 


Approximately 4-50 days of £ 0,5 MeV-energy proton observations starting 
In early 1983 appear in Figure 9 To facilitate comparisons, the observation 
times of spacecraft data have been translated back to the sun assuming a sun- 
spacecraft propagation delay corresponding to a corotating flow with a con- 
stant speed of 500 km/s and are plotted against Carrington Solar Rotation 
number. An average radial propagation speed of 500 ± 35 km/s was derived us- 
ing spacecraft radial and longitudinal positions, a solar rotation period of 25.38 
days, and measured delay times between correlated recurrent features. Since 
the spacecraft were close in longitude (<23°), the small additional adjustment 
for corotation delay was unnecessary in the figure. Features associated with a 
CIR characterized by 500 km/s line up vertically at the three spacecraft. Re- 
current features at Pioneer 11 and Voyager 1 which are associated with a 
higher speed corotating structure appear to the left of a vertical alignment 
with Voyager 2, slower characteristic speed structures appear to the right. 
Similarly, latitude variations of the boundary and width of a solar wind stream's 
source region cause fluctuation? around this vertical organization. 

Temporal periodicity 13 apparent at each spacecraft throughout this time 
period. This is in direct contrast to observations from Pioneer 10 during the 
rise to maximum solar activity, when no recurring features were observed 
beyond 14 AU. During the first twelve solar rotations shown In the figure, the 
average periods of the recurring proton flux enhancements in the spacecraft 
reference frames are ~28 days. 

Double peaks on a number of recurrent flux enhancements may indicate 
the presence of either the forward- reverse shock pair on a single CIR or single 
shocks on two or more merging CIRs. 

Three fairly distinct intervals are apparent in Figure 9, During Interval 1, 
lasting until mid-rotation 1737, the counting rates at all three spacecraft are 
highly correlated and are dominated by a common 26-day recurrence. The re- 
current nature of the proton flux enhancements indicates that CIRs are most 
likely present, enforcing an ordered structure throughout this region of the 
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Fig 9: Upper panel: counting rates of £ 0.5 MeV-energy protons from tae Voy- 
ager 1 and 2 and Pioneer 11 spacecraft are plotted versus observation time 
(Voyager 2 time is indicated) translated back to the sun and identified by Car- 
rington solar rotation number Dashed lines separate the three analysis inter- 
vals discussed in the text. 

Lower panel: He 10830 X coronal hole boundaries are plotted versus Carrington 
solar rotation number 

heliosphere The salient point of these observations is the periodicity of the re- 
currence, not an exact correspondence of the intensity- timo profiles. 

Interval 2, from mid-rotation 1737 to mid-rotation 1744. is identified by 
counting rates at all three spacecraft lower than those observed during Interval 
1. A common 26-day recurrence is still apparent, although, at times the peaks 
are barely visible at Voyager 1. On the average, the relative decrease in peak 
counting rates at Voyager 1 and at Pioneer 1 1 is approximately twice the rela- 
tive decrease at Voyager 2. Since Pioneer 11 is closer to Voyager 1 in helio- 
graphic latitude but closer to Voyager 2 in distance and longitude, the similan- 
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ty between Pioneer 11 and Voyager 1 observations Indicates that the overall 
variation In the fluxes is most likely latitudinal, and not radial or temporal in 
nature. 

Interval 3 (after mid-rotation 1744) is characterized by peak flux levels at 
each spacecraft comparable to those in Interval 1, although now the intensity- 
time profiles are distinctly different. The interleaving of two 28-day re- 
currences results in the apparent 13-day periodicity continuing from late in In- 
terval 2 throughout Interveil 3, TTie interleaved 28-day recurrences most prob- 
ably result from distinct solar sources well separated in hellographic longitude. 
Strong 13-day recurrences appearing at Voyager 2 and at Pioneer 11 during ro- 
tations 1742 and 1743 do not appear at Voyager 1. Thirteen-day peaks during 
rotations 1744 and 1745 are strongest at Voyager 1, barely discernible at 
Pioneer 11, and absent at Voyager 2. Since temporal shifts due to radial propa- 
gation have already been removed, the staggered appearance of these 13-day 
series is another Indication that distinct latitudinal differences exist and per- 
sist into the outer solar system 

Correlations with Coronal HoIb Dynamics: ’Tgure 9 also displays estimates of 
He 10830 Acoronal ho^ boundaries plotted versus central meridian observation 
time (identified by Carrington solar rotation number). These inferred coronal 
hole boundaries and photospheric magnetic polarity information, both derived 
from earth-based observations, are Included in the Ha solar synoptic charts 
published monthly in Solar- Geophysical Data Reports. Coronal holes appearing 
in regions of dominantly north (south) solar magnetic polarity are shaded 
(unshaded). 

The dominant near equatorial coronal hole during Interval 1 is a long-lived 
(~7 solar rotations) equatorial extension of the north polar hole. Numerous 
mid-latitude and equatorial coronal holes, as well as an equatorial extension of 
the south polar hole, are also visible. Equatorial extensions of polar holes are 
most often associated with the highest speed solar wind near the heiiographic 
equator, and higher speed streams are known to dominate control of the inter- 
planetary medium. Since two of the spacecraft are north of the equator, the 
equatorial extension of the north polar coronal hole seems to be the most prob- 
able source of the corotating plasma feature controlling the interplanetary 
medium in the region of tile spacecraft. For part of Interval 2 (rotations 1736 
to 1744), the north polar extension recedes to well above 45° while two south 
polar extensions grow, reaching toward the equator, coincident with decreased 
flux levels at Voyager 1 and Pioneor 11 and with the continuing presence of 
moderate flux levels at Voyager 2. This suggests that the high speed streams at 
Voyager 1, at ,$ 20°N latitude, may emanate from mid-northern latitudes, while 
those responsible for the continuing moderate flux levels at Voyager 2 may be 
from a combination of the new equatorial extensions of the south polar hole 
and the continuing near equatorial holes of northern polarity. Finally, from ro- 
tations 1745 to 1747 during Interval 3, the north polar coronal hole re-extends 
to the equator while, at the same time, the intensities of the recurrent flux 
enhancements reach levels comparable to those in Interval 1. 

The correlated, recurrent enhancements of £ 0,5 MeV-energy proton flux 
observed at these three spacecraft for approximately six solar rotations in ear- 
ly 1983 demonstrate that Long-lasting structuring of the interplanetary medium 
by solar wind streams and their CIRs can persist over a latitude range of ~2Q a 
out to ~ 17 AU. Major differences in the intensity of recurrent proton, flux 
enhancements then continue for over ten solar rotations from mid- 1983 into 
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1964, showing that large-scale differences can also exist over a latitude range 
of only -20°. Additionally, recurrent proton flux enhancements occurring twice 
per solar rotation are observed out to at least 19 AU, indicating that CIRs do not 
necessarily merge within this radius. 

These enhancements, which depend upon the dynamics and structure of 
the heliosphere in the outer solar system, and their apparent correlation mth 
coronal hole dynamics suggest that structuring of the interplanetary medium is 
still generally controlled by coronal hole dynamics out to ~20 AU and up to 23° 
latitude from the heliographic equator. 


1.2.3. A Heavy Isotope Spectrometer Telescope (HIST) Launched on ISLfi-3 in 
August 1970 

HIST is designed to measure the isotope abundances and energy spectra of 
solar and galactic cosmic rays for all elements from lithium to nickel (3 £ 7 , £ 
28) over an energy range from several MeV/nucleon to several hundred 
MeV/nucleon. Such measurements are of Importance to the study of the 
isotopic constitution of solar matter and of cosmic l ay sources, the study of 
nucleosynthesis, questions of solar^system origin, studies of acceleration 
processes and studies of the life history of cosmic rays In the galaxy. 

HIST was successfully launched on ISEE-3 and provided high resolution 
measurements of solar and galactic cosmic ray isotopes until December 1976, 
when a component failure reduced its isotope resolution capability. Since that 
time, the instrument has been operating as an element spectrometer for solar 
flare and interplanetary particle studies, 

Our work on solar flare, Interplanetary, and galactic cosmic ray isotopes 
has resulted tn the following recent papers. 


• "A High Resolution Study of the Isotopes of Solar Flare Nuclei, 1 ' R. A. 
Mewaldt et a!., Ap. J. 200, 892 (1984). 

• "The Isotopic Composition uf the Anomalous Low-Energy Cosmic Rays," 
R.A. Mewaldt et al„ Ap. J. 283, 450 (1984). 

• "Measuements of Fe and Ar Fragmentation Cross Sections," K. H. Lau 
et al., proc. 19th International Cosmic Ray Conference, La Jolla, California 
(1985). 

• "Solar Cycle Variations of the Anomalous Cosmic Ray Component," R. A. 
Mewaldt and E, C. Stone, Proc. 19th International Cosmic Ray Conference, 
La Jolla, California (1985). 

•"Cosmic Ray 3 He Measurements," R. A. Mewaldt, Proc. 19th International 
Cosmic Ray Conference, La Jolla, California (1985). 

A summary of some of this work appears below. 

A New Measurement of the Isotopic Composition of Fe Fragments 

One of the natural by-products of calibrating a cosmic ray spectrometer 
with heavy ion beams is a set of data that can be used to investigate nuclear 
fragmentation cross-sect.ions of interest to cosmic ray propagation studies. 
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Our 1SEE-3 Instrument was calibrated at tbe Uevalac with beams of afl Fe 
Incident on a CHp target at ~ 500 MeV/nucleon. Products of fragmentation 
reactions In the target were analyzed by HIST as part of the calibration 
procedure for determining the isotope "tracks" in the various detectors of the 
instrument. Although the experimental satup did not permit absolute cross- 
section measurements, it did permit an accurate measurement of the relative 
yield of the isotopes for a number of elements including Cr. Mn. Fe, and Co. 
This work constitutes one part of the PhD thesis of K. H. Lau. 

An example of the mass distributions obtained for Mn and Fe Isotopes is 
shown in Figure 10. To our knowledge the observed mass resolution of ~ 0.25 
amu is the best so far achieved with a cosmic ray Instrument for Isotopes as 
heavy as a8 Fe. Figure 11 shows the relative yields of Cr, Mn, Fe and Co Isotopes 
that were observed, and compares these observations with the expected yields. 
The expected yields were circulated using a Monte Carlo calculation based on 
the semi-emplncal cross-section formulae of Tsao and SLlberberg. This Monte 
Carlo approach wr,s developed by Mark Wiedenbeck of the University of Chicago 
in an earlier coL'aDoratlve study with us of 40 Ar fragmentation, It is hoped that 
some of the differences that can be seen between the measured and calculated 
values in Figure 1 1 will serve as a basis for improvements to the semi-empincal 
cross-section formulae. 
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1.2.4. A Heavy Nuclei Experiment (HNK) Launched on HKAO-C in September 
)'J79 

The Heavy Nuclei Experiment Is a Joint experiment Involving this group and 
M. H. Israel, J, Klarmann, W, R. Binns (Washington University) and C, J. 
Waddlngton (University of Minnesota). HNE is designed to measure the 
elemental abundances of relativistic high-Z cosmic ray nuclei (17 ^ Z ^ 130). 
The results of such measurements are of significance to the studies of 
nucleosynthesis and stellar structures, the existence of extreme transuranic 
nuclei, the ongtn of cosmic rays, and the physical properties of the mterstellar 
medium. HNE vvas successfully launched on HEAO-3 and operated until late May 
1981. 

The following talks and papers were presented during the reporting period: 

• "Elemental Abundances of Ultraheavy Cosmic Rays," W.R. BLnns et ai„ 
Adv. Space Res. 4, 25 (1984). 

■ "HEAO-3 Observations on Nuclei of Platinum and Lead in the Cosmic 
Radiation," D.J. Flxsen ot al., Bull, of Amur Phys. Soc. 29, 735 (1984), 

• "Rigidity Spectra of Ultraheavy Cosmic Ray Results from HEAO-3," M.D 
Jones et al.. Bull. Am. Riys. Soc. 29, 735 (1984). 

• "Elemental Abundances for Cosmic Rays with 26 £ Z ^ 42: HEAO-3 
Results," W.R. Binns et al.. Bull. Am. Phys. 5bc. 29, 734 (1984). 

The Response of Ionization Chambers to Relativistic Heavy Nuclei 

Gas-fitied ionization chambers have been used on recent cosmic ray 
experiments designed to measure heavy element abundances, such as the 
Heavy Nuclei Experiment (HNE) on HEAO-3 and Ariel VI. The response of such 
chambers is expected to be proportional to the energy deposited by the 
particle traversing them. At low energies this energy deposit is simply the 
iomzatlon energy loss, while at high energies energetic knockon electrons are 
able to escape from the chamber, reducing the energy deposit. 

To first order the Ionization energy loss scales as the square of the particle 
charge Z, however at high Z this assumption breaks down. A more complete 
expression Is given by Ahlen, and predicts an energy loss rising slightly faster 
than Z 2 . Such effects are important when identifying ultraheavy elements. 

We have performed two calibrations of ion chambers at the LBL Bevalac 
using beams ranging from 2 aMn to 92 U. The first, in 1982, was done with a 
prototype of the HNE ion chamber module which was essentially identical to 
those used in flight. Thus those data, reported in the 18th ICRC, are directly 
applicable to our flight experience at the energies calibrated. The second 
calibration, in 1964, used lab chambers which vmre made of thinner and more 
uniform materials, permitting better resolution and better knowledge of the 
beam energy in each Ion chamber, at the cost of less direct relevance to the 
flight data. Figure 12 is a schematic drawing of the 1984 detector. 
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Figure 12 

Particles entering the L984 detector traversed -0.1 gem -2 of sheet 
material in the upstream window, rather than the ~1 gem -2 of alumlmir. 
honeycomb In the flight prototype; thus the energy loss In the window is much 
smaller and more uniform. Also, in the 1904 calibration the beam energy was 
measured with a magnetic spectrometer after being degraded to the calibration 
energy, rather than being calculated from, an energy loss model. 

The lab ion chambers had aluminized mylar electrodes (0,0 mgcm" ? ) rather 
than aluminum screenwire (10 mil diameter, 02.5 mil spacing); thus the 
production of knockons la much more uniform. A Monte Carlo model of 
knockon production correctly predicts the degradation in resolution caused by 
non-uniform production of knockons tn the screen wire electrodes. This 
resolution degradation in the flight chambers tends to mask the relatively 
subtle deviations from Z 2 scaling. 

In 1982 the ultraheavy capabilities of the Bevaiac were new and we 
calibrated only on beams of -1700 MeV/arnu 2s Mn and ~1000 MeV/amu 79 Au. in 
1904 beams of 20 Fe, 38 Kr. a*Xe, 67 Ho, 7fl Au. and 92 U nuclei were available at 
maximum energies of 1606 MeV/amu for Fe to -039 MeV/amu for U. 



Fiaure 13 
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Result* of the 1964 Calibration. Flguro 13 ghowa the response of chambers 1,5, 
and 0 lo nuclol aa a function of the energy at the mldplano of the 

appropriate chamber, and compares thoLr signals to Lho calculated dE/dx, 
arbitrarily nomializod at 500 MoV/amu. This normalization corresponds Lo 27 9 
eV per Ion pair at low energies for the P-10 gas used (90% argon. 10% methane) 

It Is apparent that the signals fall below that predicted by dE/dx at 
energies above 700 MoV/amu, This Is somewhat surprising since at those 
onorgles wo would expect knockons escaping from the exit window Lo bo in 
equilibrium with those arriving from above, particularly for chambers 5 and 0 
which have ~2 gem’ 2 of upstream material. However, wo note that a knockon of 
a particular energy will be scattered closer to 90" from the incident particle 
direction as the Incident energy Increases, and that for energetic knockons the 
scattering length in a gas may be several cm. Thus the decrease in ohservod 
signal may be partly due to knockons escaping from the sides of the chambers, 
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By Interpolating to a particular energy wo can construct a plot of signal 
vorsus Z at that energy At low energies, the heaviest nuclei have an effective 
charge, Z = Z[ 1 -expf-lOO/JZ -273 )]. due to electron capture. Figure 14 shows 
the pulao heights, scalod down by ZjS«. at four energies for Z = 26-92, using ion 
chambers 1,6, and 6. The straight linos represent a linear fit to the data 
excluding U, and It Is apparent that there Is a small nogatlvo non-Z 2 effect. The 
chargo of an flZ Pb nucleus would be underestimated by about 0 6 chorge units 
at those energtos, in contrast with the chargo overestimate of t-3 chargo units 
obsorvod for the flight chambers. 

Although the non-Z z effects Ln these chambers differ from Lhoso obsorvod 
tn the prototype flight chambers, the assumption of Z 2 scaling is still not 
seriously in error. Wo also note that our published abundances above chargo 50 
have used the Ceronkov detector alone to assign charges, and are unnffectod by 
small non-Z B effects In the ion chambers. 

Since the two calibrations differ, the Ionization response to energy loss 
must be sensitive to dotails of the mass distribution above, below, and within 
the chambers. As a result we are using the flight data to directly determine 
both the energy dependence and effective non -7r correction. 


1 .2.6. Models of Saturn's Magnetic Field 

A reanaiysia of the Vector Kalium Magnetometer data taken by Pioneer 11 
during its encounter '"ith Saturn Ln 1979 has provided a very good 
axisymmetrtc model (designated as the P u 84 model) of the magnetic field 
within 6 Saturn radii of the planet. The appropriately weighted root mean 
square average of the difference between the observed and the modeled field is 
1,13% for all the data Inside 0 Saturn radii and is 0.40% for all data inside 3 
Saturn radii, there being no data inside 1.35 Saturn radii. The difference inside 
8 Saturn radii Is 1.81% If the P u 64 model Is replaced by the Z 3 model of 
Connerney, Ness, and Acuna, which Is based on data taken by Voyagers 1 and 2. 
The difference Is 1.02% If the Z 3 model Ls used for the Pioneer 11 daLa inside 3 
Saturn radii. Both the P;jfH and the Zg models are good models, but the P n 84 
model Is somewhat better. 

The P n 84 model Ls better than previous models based on the Pioneer 11 
data for several reasons. L. R. Doose has recently provided more accurate data 
on the roll attitude of Pioneer 11 during the Saturn encounter based on a 
reexamination of the Pioneer 11 Imaging Photopolarlmetar data of GehreLs ot 
al. Some of the one minute averages of the Pioneer II data had been excluded 
from previous analyses because of concerns of then- reliability but a careful 
review determined that they were of high reliability and they are now included 
in the analysis. It was found ti-at the currents In the bow shock, magnetopause, 
and ring current regions varied during the encounter and hence that the 
coefficients in the model that specify the field due to these external sources 
must have different values after the data gap during occultatlon than they had 
before. 

For the P n 84 model the internal source coefficients are g? = .2114. g§ = 
.0180, gS - .0226 Gauss for the internal sources; and for external source 
coefficients G? = -7,lnT while Pioneer 11 was approaching Saturn and was -12.6 
nT while It was leaving. 
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1.2.6. Proposal for an Advanced Conpodtion Explorer (ACE) 

This Investigation, proposed Jointly by this group, and by W. D. Arnett and J. 
A. Simpson (University of Chicago), L F. Burlaga (GSFC). R. E. Gold and S. M. 
Krimigts (APL/JHU), W, C. Feldman (LANL), G. Gloookler and G. M. Mason (UMd), 
and J. V Hollweg (UNH), Is for the study of an Advanced Composition Explorer 
(ACE). This Explorer-class mission would moke comprehensive measurements 
of the elemental and Isotopic composition of accelerated nuclei with Increased 
sensitivity of several ordors of magnitude, and with tmprovod mass and charge 
resolution. ACE would observe particles of soLar, interplanetary, and galactic 
origins, spanning the energy range from that of the solar wind (~ 1 
keV/nucleon) to galactic cosmic ray energies (several hundred MeV/nucleon). 
Definitive studios would be made of the abundance of essentially all Isotopes 
from H to Zn (lsZs:30), with exploratory isotope studies extending to Zr (Z=40), 
and element studios extending to U (Z=92), 

ACE would bo a coordinated experimental and theoretical effort, designed 
to investigate a wide range of fundamental problems. In particular, ACE would 
provide the first extensive tabulation of solar isotopic abundances based on 
direct sampling of solar matorlal and would establish the pattern of isotopic 
differtmcas between galactic cosmic ray and solar system matter Those 
composition data would be used to investigate basic dynamical processes that 
Include the formation of the solar corona, the acceleration of the solar wind, 
and the acceleration and propagation of energetic nuclei on the Sun, In 
Interplanetary space, and In cosmic ray sources. They would also be used to 
study the history of solar system material and of galactic cosmic ray material, 
and to Investigate the differences In their origin and evolution. 

The ACE study payload Includes four high resolution spectrometers, each 
designed to provide the ultimate charge and mass resolution In Its particular 
energy range, and each having a collecting power l to 3 orders of magnitude 
greater than previous or planned experiments. Included In the study would be 
two spectrometers, a SoLar Isotope Spectrometer (SIS) and a Cosmic Ray 
Isotope Spectrometer (CRIS), for which Caltech would play a leading role. These 
spectrometers would make use of the proven mass-resolution techniques and 
Large-area detectors that were developed and tested by this laboratory over the 
past decade, partly through the support of this grant. 


1.2.7. Galileo Heavy Ion Counter 

This experiment, being constructed by this group in collaboration with N. 
Gehrels at Goddard Space Flight Center, Is being added to the Galileo mission as 
an engineering subsystem. It will monitor penetrating (~ ID to ~ 200 
MeV/nucleon) sulfur, oxygen, and other heavy elements In the Jovian 
magnetosphere with the sensitivity needed to worn of potential "single-event 
upsets" (SEU) In the attitude control system computer. (SEUs are state 
changes induced by Ionizing radiation.) Caltech i3 responsible for management, 
detector testing, and calibration of the experiment, which Is based on 
repackaging the Voyager CRS prototype unit (the PTM). Although the primary 
purpose Is engineering support, the data will allow us to continue our 
Investigation of spectra of trapped Ions in the Jovian magnetosphere and their 
relation to the Jovian aurora. In addition, during cruise phase and in the outer 
Jovian magnetosphere, we will use the Instrument to measure the elemental 


composition of solar Haro ovonts and of the anomalous cosmic ray component, 
beginning In 1906 and continuing through ~ 1991. 


2. Gamma Ray? 

This rosoarch program, which has roculved significant support from 
Caltech, is directed toward the Investigation of galactic, oxtragalactlc, and solar 
gamma rays with spectrometers of high angular resolution and modorato 
energy resolution carriud on spacecraft and balloons. The main oflorta havo 
been directed toward the following two categories of experiments. 


2.1. Activities in Support of or in Preparation for Spacecraft 
Experiments 

These activities generally embrace prototypes of experiments on existing 
or future NASA spacecraft and they complement and/or support such 
experiments. 


2.1.1. A Balloon- Borne Gamma Ray Imaging Payload (CROP) 

The principal focus of the current gamma- ray astronomy effort is the 
completion of a balloon-borne imaging gamma-ray telescope for galactic and 
oxtragalactlc astronomy observations. A shielded Nal Anger camera is used in 
combination with a 2-cm thick, lead, rotating, coded-aperture mask to achieve 
an imaging capability of 1070 0 6° pixels in a 20° field of view and a localization 
capdbilily of 3 arc minutes for 10c sources. This performance represents more 
than an order of magnitude Improvement over previous balloon and satellite 
instrumentation. The 16" x 2" Nal Anger camera plate will have an energy 
range of 30 keV to 5 MeV and achieve a continuum sensitivity of 
0x 10~ 7 photons/ cm 2 skeV at 1 MeV 

Figure 15 s_ows a drawing of the detector and coded aperture mask 
systems of the GRIP Instrument, The table lists the key parameters of the 
instrument. Figures 16a and b are photographs of the instrument taken in 
early April. 




Figure 16a 


Figure 16b 
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GRIP Balloon Telescope 

Primary Detector 

40 cm x 5 cm NaJ Anger Camera 

Position Resolution <5 mm FWHM above 10C keV 

Shield 

Back Plate: 5 cm Nal 

Side 18 cm plastic scintillator 

Mask 

Hexagonal URA: 1.1m diameter x 2 cm Pb 
Rotation rate: 1 rpm 
Spacing 2.5 m from Nal detector 
No of Cells: 1750 (2.5 cm) 

j Energy Range 

0 03 - 5 MeV 

Energy Resolution 

8 3 keV O 50 keV. 70 keV 9 1 MeV 

Sensitivity 
(3a - 8 hr) 
(Southern 
Hemisphere) 

Continuum: (AE/E = l) 

9 100 keV - 1x10*° ph/cm 2 s keV 
O 1 MeV - lxl0‘ a ph/cm 2 s keV 
Broad line (AE/E = 12 x FWHM resolution) 
-3X1C 4 Dh/cra* s: 100 keV - l MeV 

i Imaging 

i 

Resolution: 1070 0.8* diameter pixels (20* FOV) 
Angular localization 3 arc min (lOo source) 


Instrument Description 

The telescope consists of a shielded detector system separated by 2.5 m 
from a lead coded- aperture mask The primary detector is a position-sensitive 
scintillator which records the characteristic spatial pattern of photons cast by 
a 7-ray source through the mask 

The mask is made of lead hexagons 2 cm thick and 2.5 cm across (flat-to- 
flat). supported by an A 1 honeycomb sandwich which is transparent at 7-ray 
energies. The aperture contains 2000 hexagonal cells of which half are open 
and half contain a lead hexagon. The cell pattern (see Figure 17) forms a 
hexagonal uniformly redundant array (HURA) that is optimal for coded- 
aperture imaging HURA's are discussed in more detail in section 2.1 2. 

Continuous mask rotation imposes an additional level of coding on the 7- 
ray signal Due to the antisymmetry of the coded-aperture pattern under 80 
degree rotation (open and closed cells interchange for all but the central cell) 
the 7-ray signal at each position on the detector is time- modulated with a 50% 
duty cycle. This feature allows a complete background subtraction to be 
performed for each detector position once every 20 seconds assuming a 1 rpm 
rotation rate. In addition, the continuous rotation permits extension of the 
field of view to 20 degrees, greatly increasing the number of pixels imaged. 

The primary detector is a Nal(Tl) camera plate 41 cm in diameter and 5 cm 
thick. The Nal is viewed by nineteen 3 inch photomultiplier tubes (PMT's) which 
are individually pulse height analyzed. The PMT gains are calibrated 
continuously using pulsed LED's for short term relative gain calibration and an 
24f Am tagged 7*ray source for long term absolute gam calibration. The Z41 Am 
source is situated 1 m above the coded- aperture mask and can be imaged 
continuously during flight, allowing a thorough checkout of the mask-detector 
imaging system 
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Background suppression is provided by an anti-coincidence shield On the 
side are 12 plastic scintillator modules which form a cylinder -10cm thick 
Each module is viewed by a single 5 inch PMT The lower shield section is a Nal 
camera plate identical to the primary sensor Further background suppression 
is provided by the primary scintillation camera itself The PMT pulse heights 
contain information on the depth of the interaction in the crystal Thus the 
lower hall of the detector can be used as an effective "lltegral shield'' for the 
reduction of background at low energies 

The telescope is mounted on an elevation pointing platform suspended 
from an azimuthal torquing system Azimuthal stabilization ar.d orientation are 
achieved using active magnetometer feedback to the azimuthal torque motor 
Elevation orientation is under command control Two magnetometers provide 
2-axis aspect information accurate to 1 5 ar " minutes This aspect information 
is recorded in the telemetry stream and allows correction of the event positions 
for pointing inaccuracies such as displacement and rotation of the telescope 
field of view 

For the initial flight of GRIP, we will record all nineteen 12 bit PMT pulse 
heights for each event Event rates of up to hxlC^/s are possible and 
consequently a data recording system with a 1 Mbit/s data rate is required We 
have developed a 14 Mb's recording system with 25 Gbyte capacity using 
commercial VCR s and audio digitizers 



Figure 17 
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The GRIP scintillation camera has been designed to have ~lcm FWHM or 
bettor position resolution over an energy range from 100 keV to 5 MeV. Figures 
10a and 10b show histograms of 7-ray event positions computed by a maximum 
likelihood method for beams of photons of 122 keV ( 07 Co) and 002 keV ( 137 Cs) 
Incident on the center of the detector. At the lower photon energy, the 10.5 mm 
FWHM of the distribution la dominated by photou statistics, At the higher 
photon energy, both Compton scattering and photon statistics contribute to 
the 7.0 mm width of the distribution, The effect of Compton scattering is most 
noticeable in the extended non-Gaussian tails of the distribution. The energy 
resolution of the GRIP scintillation camera is comparable to that of single PMT 
Nal detectors. We have measured a resolution of 7% FWHM at 562 keV. 
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Figure 18a Figure 18b 

Figures 19a and 19b are laboratory images of 122 keV and 662 keV 7-ray 
sources taken with the fully configured GRIP telescope. The sources were 
suspended 10 m from the coded-aperture mask and the imaging algorithms 
were adjusted to account for the finite distance to the source. The Images 
demonstrate the ability of the telescope to locate and resolve pomt sources at 
7-ray energies. 
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Figure 19a 


Figure 19b 


The instrument technique has been described m. 

• "Gamma-Ray Imaging with a Rotating Hexagonal Uniformly Redundant Array 
W R Cook et al IEEE Trans Ituclear Sex NS-31. 77". -77b ( 1984) 

• A Thick Anger Camera for Gamma- Ray Astronomy W R Cock et al IEEE 
Trans on Nuclear Science NS-31. 348 (1984) 

The GRIP telescope is scheduled for an initial flight from Palestine, Texas in 
Pall 1985 Observing targets include the Cygnus region. NGC4151. and the Crab 
region Future flights are anticipated from both the northern and southern 
hemispheres The status of the various instrument subsystems i5 as follows 

• Telescope Nal Anger cameras integrated with shield system, telescope 
enclosure, and flight electronics High speed VCR recording system 
taking data at 14 Mb/s Image of laboratory radioactive sources 
recorded and processed 

• Pointing Systems (Supported by Caltech Funds) Pointing platform 
fabricated Azimuthal torquer fabricated and tested Elevation motor, 
ballscrew. and encoder procured and installed 

• Support Systems Precision magnetometers procured and tested 
FTecision rotation stage procured and tested Integration of 
magnetometers and rotation stage in progress VCR recording system 
with 8 VCR units and 25 Gbyte capacity under fabrication Power system 
designed and under fabrication Temperature probe system designed 
and under fabrication 


We discuss below the coded aperture mask technique used in GRIP 
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2.1.2. Hexagonal Uniformly Redundant Armya for Coded Aperture Imag in g 

Introduction to URA'a. Uniformly redundant arrays are used In coded-aperture 
Imag ing, a technique tor forming Images without mirrors or leases, This 
technique Is especially Important for the high energy x-ray and v-ray region 
above 20 keV. In this technique, a mask consisting of opaque (closed) and 
transparent (open) areas is placed between the photon sources to be imaged 
and a position sensitive detector or a detector array, Each source casts a 
shadow pattern of the mask or aperture onto the detector. This shadow pattern 
may be viewed as an encoded signal for that source direction. If each possible 
source code is unique, the detected composite of overlapping shadow patterns 
may be decoded to produce an image of the source distribution. 

Figure 17 shows a mask suitable for imaging. This mask consists of an array 
of open (white) and closed (gray) ceils arranged in a periodic pattern. The unit 
pattern is outlined. The mask in figure 17 is a uniformly redundant array 
(URA). URAs have an especially desirable property: the overlap between two 
source codes is independent of the source directions as long as the sources are 
sufficiently separated. Except for periodicity, this guarantees a unique decoding 
of the composite shadow pattern with a maximal immunity to statistical, noise. 
To date, most work on URAs has concentrated on those constructed on 
rectangular lattices. Here we focus on URAs constructed on hexagonal lattices, 
although many of the results are independent of the lattice type. 

We have developed procedures for the Construction of a special class of 
URAs, the skew- Hadamard URAs, which have the following properties: 

1) They are nearly half open and half closed. 

2) They are antisymmetric (exchanging open and closed cells) upon rotation, 

by 180° except for the central cell and its repetitions. 

Some of the skew-Hadamard URAs constructed on a hexagonal lattice have 
additional symmetries. These special URAs that have a hexagonal unit pattern, 
and are antisymmetric upon rotation by 60°, we call hexagonal uniformly 
redundant arrays (HURAs). The mask in figure 17 is an HURA. 

HURAs are particularly suited to our application, 7-ray imaging in high 
background situations. In a high background situation the best sensitivity is 
obtained with a half open and half closed mask. Furthermore, systematic 
variations of the detector background from position to position can be larger 
than the variations in detected flux due to sources. With a skew-Hadatnard URA 
a simple rotation turns the mask into a near anti-mask, allowing exact 
position-by-position background subtraction. Also, the hexagonal symmetry of 
an HURA is more appropriate for a round position-sensitive detector or a 
close-packed array of detectors than a rectangular symmetry. This is especially 
true for shielded detector systems where compactness is at a premium 

Mathematical Structure of URAs 

A URA is defined within a unit pattern which is repeated periodically. The 
number of cells in this unit pattern is the order v of the URA. Of these cells, k 
of them are closed and v— k are open. The uniform redundancy property of 
URAs involves how frequently a given displacement between closed ceLla occurs. 
We will consider a cell within a repetition of the unit pattern as equivalent to 
the corresponding cell in the unit pattern, and will therefore define the 
difference between two cell centers as the vector displacement between them 
folded back into the unit pattern. For a URA, all possible differences occur a 
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uniform number X limes among the pairs of closed cell centers This property 
guarantees the uniform overlap of source codes discussed in the introduction 

The mathematical structure of a URA is that of an Abelian group difference 
set. which is specified by an Abelian (additive) group G of order v, and a set D of 
k elements of G with the property that any possible nonzero difference occurs 
exactly X times between elements of D For a URA the group G is the lattice 
translations modulo the periods of the mask pattern, and the set L) contains 
those translations that take the central cell to a closed cell The simplest 
examples of group difference sets are one-dimensional sets known as cyclic 
(ixffmmca sats defined on the group of integers mod v These play an important 
role in the construction of many URAs URAs in the class consided here, the 
skrw- Hadamard. URAs are nearly antisymmetric That is. for any nonzero 
element in the group G. either it or its negative but not both are contained in 
the difference set D These skew- Hadamard URAs are a subset of the Hadamard 
URAs which are nearly half open and half closed Hadamard LRAs are 
characterized by the parameters v=4n-l, k~2n-l. X = n-1 for some integer n 
Two interesting facts about skew- Hadamard URAs have been proven 

'. ) All skew- Hadamard URAs have a cyclic group G. and therefore can be 
constructed fr"*m skew- Hadamard cyclic difference sets 
2) All skew- Hadamard cyclic difference sets are of prime order v = 3 mod 4 
and can be geneiated from '.he quadratic residues mod v 
These facts allow us to present a construction for all antisymmetric or 3 kew- 
Hadamard URAs 

Construction o' Skew- Hadamard URAs from Quadratic Residues 



Figure 20 
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A sLmple procedure for generating any skew-Hadamard URA has been 
developed. An example constructed on a hexagonal lattice Is shown In figure 
20. The procedure consists of the following steps; 

1) Choose the lattice on which the URA is to be constructed. The lattice i3 
defined by picking two basis vectors 3 0 and Si . For our example we have 
chosen a hexagonal lattice, which has the basis vectors separated by 60" 

2) Choose as the order of the URA a prime of the form v = 4n-l . In our 
example v = 23 . 

3) Construct the order v skew-Hadamard cyclic difference set from the 
formula 

D=|i 2 .2 2 (—i-) 2 } mod v (1) 

4) Choose an integer r and label ail the cells so that the cell centered at 
i0 o +JSi is labeled with 

l “ (l + rj) mod v (2) 

and make all cells with labels in D closec. In our example r=5. 

The heart of this procedure is the construction of the skew-Hadamard 
cyclic difference set in step 3. Step 4 transfers the difference set properties 
onto the lattice, This is done through the labeling, which transforms addition 
mod v to vector addition on the latttce modulo the resulting lattice periods. 

The freedom available in this procedure rests in the choice of the lattice, 
the choice of the order v.and the choice of the multiplier r. The lattice type will 
determine what symmetries can occur. The possible orders form a rather dense 
set, the first few choices being v=3. 7, 1 1,19,23,31 .43,47,59,67,71.79, and 83. The 
multiplier r determines the periods of the URA, and hence the shape of the unit 
pattern. Many of the v available choices result in URAs that are related by the 
symmetries of the lattice. 

Hexagonal Uniformly Redundant Arrays 

Of the large number of skew-Hadamard URAs, ail of which can be 
constructed by the procedure in the section above, we wish to pick out those 
that have a hexagonal unit pattern when constructed on a hexagonal lattice. 
These we call hexagonal -uniformly redundant arrays (HURAs). For an HURA 
each period when rotated by 60° is again a period. It can be shown from 
equation (2) that this is only possible if the multiplier r satisfies 

r 2 = r-1 mod v (3) 

This property has a simple geometric interpretation: a cell labeled l when 
rotated by 60" will have the label rl mod v. This feature, and the properties of 
quadratic residues modulo a prime, causes this restricted srt of URAs to have a 
rotational antisymmetry upon rotation by 60° as well as 180°, 

In figures 21a thru 21c we show examples of HURAs with orders 139, 151, 
and 331. 
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2.2. Experiments on NASA Spacecraft 

The scientific vitality of the SRL gamma ray program ts enhanced by 
collaboration In the analysis of gamma ray data from recent missions. This 
analysis effort is carried out on the SRL data analysis system. 


2.2.1. The Gamma Ray Spectrometer Experiment on the Solar Maxi rm im Mission 
(SUM) 

The Solar Maximum Mission satellite was launched in early 1900 and has 
been In continuous operation observing gamma-ray emission from solar flares. 
The University of New Hampshire Instrument la sensitive to photons In the 
energy range 0.3-100 MeV, and can also detect energetic neutrons above 40 
MeV. 

We have collaborated with U.N.H. Ln the study of 2.22 MeV emission from 
neutron capture on hydrogen, following up on our earlier work with the HEAO-3 
spacecraft. The focus of the SMM data analysis has been a study of 2.22 MeV 
line emission In a sample of 8 flares, which Includes the very Intense solar flare 
of 3 June 1902. The data from these flares have allowed a detailed study of the 
time history of the 2.22 MeV neutron capture line from which conclusions have 
been drawn concerning the production of low energy neutrons In solar flares, 
the density at which neutrons are captured, and the 3 He abundance In the 
photosphere. 

The following invited talk was presented in connection with this work: 

• “Gamma-Ray and Neutron Emission from Flares," T. A. Prince, Bull. Am. 

Astro 5bc 16. 51 i (1904). 

• "The Energetic Particle and Photon Component of Solar Flares." T.A. 

Prince, Bull. Am. Phys. Sac. 29, 71 (1904). 


2.2.2. Development of a /-Ray imager with Arc Second Resolution 

A major new activity In the last year has been the definition of an ultra- 
hlgh-resolutlon 7-ray Imager for solar and cosmic 7-ray astronomy. Above 100 
koV, past Instruments have been limited to a typical angular resolution of 10° 
or poorer. The GRIP project described above Is a significant step forward with 
an angular resolution of 0.5° up to energies greater than 1 MeV. It became 
clear during the last year that an extension of the GRIP technology allows 
imaging with a resolution of 1 second of arc. A 7-ray imaging device with 1 arc 
second resolution would be a 10,000 fold improvement over conventional non- 
imaging Instrumentation and have substantial new capabilities for observations 
of astrophysical gamma-ray sources. With a resolution comparable to that of 
the VIA and the Einstein X-ray Observatory, the 7-ray Imager would make 
significant contributions to our knowledge of solar flare physics, the galactic 
center, AGN, and other hard x-ray/7-ray sources. 

The arc second 7-ray imager Is based on two previously developed 
concepts. Fourier transform techniques proposed for x-ray imaging on NASA's 
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Plnhole/Occultor Facility (P/OF) can be combined with the position sonsitlvo 
dotactor development carried out by our group for the GRIP bailoan project to 
achieve greatly tmprovod imaging performance. In ossence, pairs of widely 
separated collimator grids are used to form a largo scale modulation pattern of 
the 7 -rays that can be measured by 7 -ray detectors with moderate spatial 
resolution (~1 cm FWHM) such as the GRIP Anger cameras. Each pair of 
collimator grids moasures the phaso and amplitude of a particular Fourier 
component of the source angular distribution. In analogy to radio astronomy, 
each pair of collimator grids provides one "baselino" for measurement of the 
angular source distribution. 

Planning for the practical Implementation of the hlgh-rcsolutlon 7 -ray 
Imager has taken place within the high-energy solar physics community This 
was a natural approach since the concept evolved in part from the P/OF project 
whose primary objectives aro hlgh-resolutlon studies of the sun. Initial 
discussions took place in October, 1904 at tho Workshop for a High-Energy 
Facility for SoLar Physics. Concurrently, the possibility of an SMM Instrument 
changoout was being discussed for the next solar maximum High resolution 7 - 
ray Imaging Is a very attractive optLon for such a mission. We investigated this 
possibility in detail and the results wore presented at two meetings at the 
Goddard Space Flight Center in February and April, I960. The result was a 
concept for a MAX ’91 mission which would include an instrument called GRID 
(Gamma-Ray Imaging Device) Caltech played a central role In the development 
of this instrument concept in close cooperation with Marshal Space Flight 
Center, Goddard Space Flight Center, UC3D, UC Berkeley, and European 
collaborators from Birmingham, U.K. and Utrecht, Netherlands. 

The GRID Instrument consists of two collimator grid planes viowed by a 
complement of x-ray and 7 -ray detectors. The collimator grid planes are 
separated by 6.7 meters. The finest collimator grid scale is 50fj. which yields an 
angular resolution of 1.5 arc seconds. Fabrication of the 1-2 cm thick 
collimators will be accomplished through stacking of commercially photoetched 
tungsten grids. For x-ray detectors, xenon-filled multiwire proportional 
chambers will be used similar to the Birmingham University Spacelab II design, 
The 7 -ray detectors will be similar In design to the Anger cameras developed for 
the GRIP balloon Instrument, The design of the GRID Imager Is modular with a 
possibility of accommodating 16-64 collimator pairs. 

In the following months, we expect to play a central role in a detailed 
Phase A study for the MAX '91 mission. Past study activities have received 
support from this grant as well as from the NSF Presidential Young Investigator 
Program 


3. Other Activltiea 

R. A. Mewaldt is serving as a member of the High Energy Astrophysics 
Management Operations Working Group (HEAMOWG), the Cosmic Ray Program 
Working Group (CRPWG), and as the "OG" program chairman for the 19th 
International Cosmic Ray Conference. 

T. A. Prince has received a Presidential Young Investigator Award from the 
National Science Foundation. 

E. C. Stone continues to serve as NASA’s Project Scientist for the Voyager 
Mission. He is also a member of the Space Science Board, the NASA University 
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Relations Study Group, and the stcoring group of the SS0 study-, "Major 
Directions for Space Sclonco: 1996-2015" 
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